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ABSTRACT 

Inspirals and mergers of black hole (BHs) and/or neutron star (NSs) binaries are expected to be 
abundant sources for ground-based gravitational-wave (GW) detectors. We assess the capabilities of 
Advanced LIGO and Virgo to measure component masses using inspiral waveform models including 
spin-precession effects using a large ensemble of GW sources randomly oriented and distributed 
uniformly in volume. For 1000 sources this yields signal-to-noise ratios between 7 and 
200. We make quantitative predictions for how well LIGO and Virgo will distinguish between BHs 
and NSs and appraise the prospect of using LIGO/Virgo observations to definitively confirm, or reject, 
the existence of a putative “mass gap” between NSs (m < 3 Mq) and BHs (m > 5 Mq). We find 
sources with the smaller mass component satisfying m 2 ^1-5 Mq to be unambiguously identified as 
containing at least one NS, while systems with m 2 ^ 6 Mq will be confirmed binary BHs. Binary BHs 
with m 2 < 5 Mq (i.e., in the gap) cannot generically be distinguished from NSBH binaries. High-mass 
NSs (2 < m < 3 Mq) are often consistent with low-mass BH (m < 5 Mq), posing a challenge for 
determining the maximum NS mass from LIGO/Virgo observations alone. Individual sources will 
seldom be measured well enough to confirm objects in the mass gap and statistical inferences drawn 
from the detected population will be strongly dependent on the underlying distribution. If nature 
happens to provide a mass distribution with the populations relatively cleanly separated in chirp mass 
space, as some population synthesis models suggest, then NSs and BHs are more easily distinguishable. 


1. INTRODUCTION 

Adv anced LIGO (lAbadie et al\ I20I4D and Advanced 
Virgo (jAcernese et 0,1.1120151 ) will be the most sensitive 
observatories in the gravitational-wave (GW) spectrum 
between 10 Hz to a few kHz. Binary systems comprised 
of compact stellar remnants, such as stellar mass black 
holes (BH) and neutron stars (NS) merge at frequencies 
between ^ 100 and ^ 1000 Hz and are the primary 
science target for the LIGO/Virgo (LV) network. 

GW observations will investigate these systems in ways 
not accessible to electromagnetic observations. Detailed 
observations of individual sources will provide unparal¬ 
leled insight i nto strong field gravity and the NS equa¬ 
tion of state (lEla nagan and Hinde rs 120081 iRead e t al\ 
2?MlLackev ct 0,1.11201 4HWade et R/.ll2ni4ILi et R/.ll2fria 

Sampson et al\ 120131) . Inferred characteristics of the 
compact binary population will feed back into compli¬ 


cated, ill-constrained phys ics of compac t obje ct forma¬ 
tion and binary evolution (|Abadie et al\ (I2010D and ref¬ 
erences therein). 

We investigate the capabilities of an Advanced LV net¬ 
work of detectors to constrain the individual component 
masses of a compact binary and thereby distinguish be¬ 
tween NSs and BHs. We identify the mass intervals for 
which neutron star black-hole binaries (NSBH), and bi¬ 
nary black hole (BBH) systems can be securely distin¬ 
guished. This is the first large-scale study to characterize 
compact binary posterior distribution functions (PDFs) 
including spin-precession effects over a broad range of 
masses and spins. 

We use our results to assess TV’s role in resolving the 
debate over the compact object “mass gap.” Observa¬ 
tions of X-ray binaries suggest a depletion in the mass 
distribution of compact remnants between the highest 
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mass NSs 2 Mq) and the lowest mass BHs (> SMq) 
(|6zel et al.ll201Cll : rFarr et al.ll2011[l . Inferring masses elec- 
tromagnetically is challenging and systematic errors may 
dominate. Including variable emission from the accre¬ 
tion flow in the analysis of the same X-ray binaries sys- 
tem atically finds lower ma sses and disfavors the mass 
gap (jKreidberg et al.ir2012D . 

The mass distribution of NSs and BHs has implications 
on plausible explanations for t h e core -collapse supernova 
mechanism. IBelczvnski et al\ (|20I2ll suggests the mass 
gap as observational evidence that supernovae develop 
rapidly (within 100 to 200 ms), while a “fille d gap” fa¬ 
vors l onger timescale explosion m echanisms. iKochanekI 
(|20I4fl and iClausen et al\ (|2015ll find that a bimodal 
mass distribution for compact remnants is a natural con¬ 
sequence of high mass red supergiants ending their lives 
as failed supernovae. 

GW observations have been suggested as a means of re¬ 
solving this controversy because the component masses 
are directly encoded in the signal. We use our study 
to forecast whether LV observations can confirm or fal¬ 
sify a mass gap between NSs and BHs. Our results are 
to be contras t ed wit h the contemporaneous paper from 
iMandel et al\ (|20I5ll . in which they assume a specific 
mass distribution (generated from population synthesis) 
and found it is possible to distinguish BHNS and BBHs 
and infer a gap after tens of detections. 

2. INFERRING PHYSICAL PARAMETERS FROM 
GRAVITATIONAL-WAVE DATA 

Inspiral waveforms are well understood from post- 
Newtonian (PN ) expansions of the binary dynamics (e.g. 
[Blanched 1200611 . PN waveforms enable template-base 
data analysis methods where many trial waveforms are 
compared to the data. The model waveforms, or tem¬ 
plates, are parameterized by the masses, spins, location, 
and orientation of the binary. 

GW signals from compact binaries will be in the LV 
sensitivity band for tens of seconds to minutes, evolv¬ 
ing through Ncycie ^ lO^-IO^ cycles before exiting the 
band or merging. Long duration signals place strict de¬ 
mands on the acceptable phase difference between the 
template waveforms and the signal. To leading or¬ 
der in the PN expansion, th e phase evolution depends 
only on the “chirp mass” isee lPeters and Mathewslll963ll 

Me = (^17712)^^^ (mi -I- m 2 ) where mi > m 2 are the 
component masses of the binary. The uncertainty in Me 
scales as (iVeyde)”^^^ (jSathvaprakash and Schut j 1200911 
and is thereby constrained to high precision while the 
component masses are completely degenerate. Higher 
order corrections introduce the mass ratio q = milmx 
and couplings between the intrinsic angular momentum 
of the component bodies 5'i^2 and the orbital angu¬ 
lar momenta of the system L breaking the degeneracy 
between mi and m 2 , though large correlations re main 
(iCutler and FlanaganI 1994 Poisson and Willlll995ll . 

iHannam et al\ (|201.1 ) investigated how well LV can dis¬ 
tinguish between NSs and BHs by approximating param¬ 
eter confidence intervals using the match between a pro¬ 
posed signal and a grid of templates. Their study used 
a simplified version of the full PN waveforms parameter¬ 
ized by a single “effective spin” and concluded that most 
LV detections will not provide unambiguous separation 


between NSs and BHs. 

The effective spin approximation overstates the 
degeneracy between mass and spin because it ig¬ 
nores p recession of the orb it al pla ne and spin align¬ 
ments. IChatziioannou et^ (1201511 revisited the topic 
of component-mass determination with a Markov chain 
Monte Garlo (MGMG) analysis using waveforms that in¬ 
cluded spin and orbital precession . Their study found 
that spin precession reduces mass-spin correlations, con¬ 
cluding that BNS systems with components consistent 
with known NSs in binaries would not be misidentified 
as either low-mass black holes or “exotic” neutron stars 
which they define as having masses either below 1 Mq or 
above 2.5 Mq, and/or with dimensionless spins x > 0.05. 

In this work we use the same methodology of 
IChatziioannou et al\ (l2015| l - employing an MCMC with 
templates that include full two-spin precession effects to 
infer the PDF - but use a large ensemble of NSBH and 
BBH signals, paying particular attention to LV’s ability 
to identify sources which have one, or both, components 
in the mass gap. 

3. METHOD 

Our study uses the LALInference software library for 
recovering the parameters of compact binary systems. 
In our work we elect to use the MCMC implementation 
lalinf erencejneme though results do not depend on the 
chosen sa mpler. A complete d escription of the software 
is found in iVeitch et al\ (1201511 . 

LALInference is part the LSC Analysis Library (LAL) 
which has a wide variety of template waveforms avail¬ 
able. For this -work we use the SpinTaylorT2 wave¬ 
form implemented in LAL. A detailed description 
of the wav e form can be found in Appendix B of 
iNitz et ali (|2013ll . To simulate a population of plausi¬ 
ble LV detections we draw 1000 binary parameter combi¬ 
nations from a uniform distribution in component masses 
with nil > W 2 , mi > 1 Mq and mi -|- m 2 < 30 Mq. The 
maximum total mass of 30 Mq is chosen so that the 
merger and ring-down portion of the waveform (where 
the PN approximation is invalid) does not dominate the 
signal as observed by LV. Dimensionless spin magnitudes 
X are drawn uniformly from [0,1] and the spin vector 
orientation is randomly distributed over a sphere with 
respect to L. Notice that our choice of simulated signals 
includes the possibi lity of neutron star s with anomalously 
high spins x > 0.7 ()Lo and Linll2011[l . 

The sky location is distributed uniformly over the ce¬ 
lestial sphere, orientations are randomly distributed, and 
the distance to the binary is uniform in volume. We 
reject sources which do not have signal-to-noise ratio 
SNR > 5 in two or more detectors. Figure [1] shows the 
mass distributions of our population. The bottom panel 
is a scatter plot of mi and m 2 colored by the source 
SNR over the Advanced LV network. It is important 
to note that our simulated population includes bi¬ 
naries at sufficiently high mass (M > 12Mq) that 
the merger and ringdown portion of the signal is 
detectable (|Buonanno et ^l2009h . Our study is 
focused on quantifying how well masses can be 
determined purely from the inspiral part of the 
waveform, for which we have reliable waveforms 
appropriate for parameter estimation (i.e., valid 
across all parameter space). Using inspiral wave- 
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forms for both our signal simulations and parame¬ 
ter recovery allows us to assess the inspiral effects 
separate from merger-ringdown effects. However 
we are in urgent need of processing waveforms 
that include the merger/ringdown and are valid 
over the full prior range to avoid systematic er¬ 
rors in analyses of real data. Currently no such 
waveforms are available. Available inspiral mod¬ 
els that include precession and merger ringdown 
must be calibrated to numerical relativity sim¬ 
ulations and to date are only va lid at mass ra¬ 
tios g > 1/4 (|Hannam et al.1 120141 1. Such limited- 
validity waveforms lead to non-quantified, sys¬ 
tematic biases and cannot be used for parameter 
estimation at present. In the absence of generic 
waveforms it may be necessary to apply a low- 
pass filter on real data to mask any merger and 
ringdown signals. 

Because GW emission is strongest along the or¬ 
bital angular momentum direction, the detected 
binaries from a population uniform in orientation 
and volume is biased against systems whose or¬ 
bital plane is edge-on to the observer (defined 
as having an inclination angle near 90°). Fig. 
shows the cumulative distribution function of our 
observed population (red) compared to an un¬ 
derlying population distributed uniformly in ori¬ 
entation. This selection effect has an impor¬ 
tant role in mass measurement for LV observa¬ 
tions. Precession-induced effects on the wave¬ 
forms, which are relied upon to break the degen¬ 
eracy between mass ratio an d spin, are less de - 
tectable for face-on systems (|Vitale et aZ.I [201411 . 
Previous studies which have used hand-selected 
populations of signals to study the effects of pre¬ 
cession have not accounted for this selection bias, 
overemphasizing the role precession can play for 
gen eric signals. For example, th e inclination used 
by IChati^ioannou et al.\ (|2015fl (63°) is larger 
than ^ 80% of our sources so the improvement 
they found will be fully realized for a small frac¬ 
tion of our simulated population. 

For each binary we com pute the re s ponse of Advanced 
LV at design sensitivity lAasi et ^ (120131) with a low 
frequency cutoff at 20 Hz. Each simulated “detection” 
is then analyzed with lalinference_mcmc which returns 
independent samples from the PDF for the model pa¬ 
rameters. We do not simulate instrument noise in this 
study because our focus is on the degree to which LV’s 
frequency-dependent sensitivity, and the flexibility of the 
PN waveforms, limit mass measurements. Adding simu¬ 
lated noise to our signals introduces uncontrollable con¬ 
tributions to the PDF without adding any value to our 
assessment of each simulated signal. Posterior distribu¬ 
tions in true detections will be altered by the particular 
noise realization in which the signal is embedded. 

4. RESULTS 

4.1. Distinguishing black holes from neutron stars 

We take m = 3 Mq to be the dividing line between 
the masses of NSs and BHs, with all NSs located below 
this threshold and all BHs located above. We will use 
5 Mq as the minimum mass of a BH when assuming the 
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Fig. 1.— The mass distribution of our simulated population 
of compact binaries in units of Mq. [Top panel] Histogram 
for mi (red, solid) and m 2 (blue dotted). The bottom panel 
is a scatter plot of mi and m 2 colored by the source’s SNR 
over the advanced LV network. 



Fig. 2.— Cumulative distribution function for detected in¬ 
clination angles (red) compared to the isotropicly oriented 
underlying population. There is a significant selection ef¬ 
fect in favor of face-on systems (inclination ~ 0°) because 
GW emission is strongest along the direction of the orbital 
angular momentum vector of the binary. Modulations in the 
waveform from spin-induced precession are less detectable for 
face-on systems. This selection effect will suppress the ability 
for precession to break the degeneracy between spin magni¬ 
tude and mass ratio. 

existence of a mass gap. The distinction between NSBHs 
and BBHs in our sample is determined by the measure¬ 
ment of the smaller mass, m 2 . Figure |3] shows the 90% 
credible intervals of the PDF for m 2 as a function of the 
true value from the simulated population. Each entry is 
colored by the SNR of the signal in the 3-detector net¬ 
work. Horizontal gray dashed lines denote the mass gap. 
We find all of our simulations with m 2 ^1-5 Mq to be 
clearly identified as containing at least one neutron star, 
however our population yields only ten such systems so 
these are small-number statistics. Most binaries with a 
NS below 2 Mq constrain the NS to have a mass be¬ 
low 3 Mq although occasional the posterior supports the 
smaller object being a low-mass BH . Detections of larger 
mass neutron stars (2 < m 2 < 3 Mq) have 90% credible 
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intervals which consistently extend into the BH regime. 
The tendency for recovery of high-mass neutron stars in 
NSBH systems to be consistent with low-mass BHs poses 
a challenge for determining the maximum NS mass from 
LV observations alone. 

At what point can we rule out the possibility that the 
system contains a neutron star, and definitively declare 
that we have detected a binary black hole? In Figure[3]we 
see that the true mass of the smaller object must exceed 

6 Mq before the 90% credible intervals rule out a NS. 
Depending on details of spin alignment and ori¬ 
entation, systems with m 2 as low as 4 Mq can be 
unambiguously identified as BBHs. Within this 
range the m 2 posteriors for BBH sources seldom reach 
below 2 Mq, so if a maximum NS mass were indepen¬ 
dently confirmed to be consistent with current observa¬ 
tions TV’s classification of NSs and BHs would improve. 

It may seem surprising in Fig [3] that the width 
of the credible intervals do not exhibit the 1/SNR 
scalin g predicted by Fisher matr ix approxima¬ 
tions (| Cutler and FlanaganI Il994fl . The Fisher 
approximation is only suitable at sufficiently high 
SNR that the posterior distribution function is 
well approximated by a multivariate Gaussian, in 
which case the inverse Fisher matrix is the co- 
variance of the posterior. Implicit in this con¬ 
dition is the assumption that the model wave¬ 
form is a linear function of the source parame¬ 
ters. See iVallisneidI ()2008n for a thorough decon¬ 
struction of Fisher matrix-based intuition being 
applied to GW signals. For typical LIGO/Virgo 
binaries these conditions are not satisfied for the 
mass parameters. The width of the m 2 credible 
intervals is driven by uncertainty in q which, due 
to the degeneracy with spin, is extremely non- 
Gaussian and can span the entire prior range. 
The chirp mass, on the other hand, is a suffi¬ 
ciently well constrained parameter at the SNRs 
in our simulated population. In Fig. we show 
the fractional 1 — a uncertainty in A4c inferred 
from the Markov chains as a function of the 
source value, with each event colored by the net¬ 
work SNR. The expected 1/SNR dependence is 
apparent for this parameter. Notice also that 
the chirp mass errors grow with increasing total 
mass. Higher mass binaries are in band for fewer 
GW cycles which directly impacts measurability 
as discussed in Sec. [2l 

4.2. Identifying systems in the mass gap 

Following iMandel l)2010f) we use our simulated LV de¬ 
tections to infer the relative fraction of NSs, BHs in the 
mass gap, and BHs above the mass gap. Because of the 
large mass-measurement uncertainties the number of de¬ 
tections needed to conclude BHs inhabit the mass gap 
is highly dependent on the underlying mass distribution. 
Depending on whether the observed population is heav¬ 
ily dominated by low- or high-mass systems, we find be¬ 
tween ten and many hundreds of detections are necessary 
to conclude (at three-sigma confidence) the gap is popu¬ 
lated. 

Our initial set of simulated signals features many high- 
mass ratio binaries. We simulated an additional 100 
sources with both objects having 3 < m < 5 Mq to check 
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Fig. 3. — 90% credible intervals for recovered m2 as a func¬ 
tion of the true mass. Each entry is colored by the net¬ 
work SNH^_oOthe_souree;_Ilon2on^ lines denote the mass 
gap (lOzel et fflZ.ll201Cll : IFarr et ^120111 1 . Credible intervals for 
high mass systems are limited from above by our prior on the 
total mass M < 30 and the m2 < mi convention. Axes are 
in units of Mq . 


0.06 

0.05 

0.04 

' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1*' 1 ' T ■ 1 ' 

_ • • • _ 

•1 

- -f/i - 

- 


1 0.03 

1 

0.02 



0.01 

• * 



, i.i *• i‘, 1 , 1 . 1*. 1^ 1 , 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 

injected 


Fig. 4. — Fractional 1 — cr error in chirp mass as a function 
of the true value. Each entry Is colored by the network SNR 
of the source. The expected SNR scaling of the errors is 
evidence in the chirp mass recovery, while it is absent in the 
m2 measurement due to the degeneracy between mass ratio 
and spin. The x-axis is in units of Mq. 

whether comparable-mass systems are easier to identify 
in the gap. The mass-gap population still suffers from 
large mass errors with > 95% of the sources having pos¬ 
terior support for a NSBH system. Figure [5] shows the 
distribution of the 90% credible interval widths for mi 
(red, solid) and m 2 (blue, dotted) of the gap sources. 
The majority of plausible mass-gap sources yield credi¬ 
ble intervals that are similar to or exceed the width of 
the mass gap. 

However, ~ 25% of the mass gap sources’ m 2 posteriors 
do not reach below 2 Mq. While we will not be able to 
say with any certainty that an individual source occupies 
the mass gap, we can often conclude that the binary 
contains either an unusually high-mass neutron star, or 
a pair of unusually low-mass black holes. 

Three of the mass-gap sources were constrained to be 
3 < m < 5 Mq. a careful investigation of these systems 
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width of 90% credible interval 


Fig. 5. — Distribution of 90% credible interval widths in Mq 
from 100 mass gap sources. The width of the mass gap is 2 
Mq, and therefore most of the BHs cannot be constrained to 
fall within the gap. 

revealed they were in low-probability alignments - all 
with the spin of the larger mass close to the orbital plane, 
and the best constrained having L nearly perpendicular 
to the line of sight. 

4.3. Restricting allowed spins for neutron stars 

All of the results thus far presented used a uniform 
prior on spin magnitude between [0,1]. The maximum 
upper limit on the sp in of a NS is y 5, 0 -7 while observed 
NS spins are lower (iLo and LinI 1201111 . We investigate 
if using a physical prior on NS spins can improve m 2 
measurement. 

We resample the posterior by imposing a maximum 
spin for component masses below 3 Mq, rejecting sam¬ 
ples from the Markov chain with X 2 > 0.7. We found 
no indication that restricting the range of neutron star 
spins significantly improves uncertainty in the determi¬ 
nation of m 2 . 

Figure |6] demonstrates the limited role of m 2 ’s spin on 
mass determination. Open circles mark the true values 
of component masses for nine representative examples. 
Going through each circle is the scatter plot of the pos¬ 
terior samples colored by X 2 ■ The arcs traced out by the 
samples are lines of constant Aic- The vertical and hor¬ 
izontal dashed lines denote the mass gap. Notice that 
there is no obvious correlation between position along 
the arc and X 2 “ the spin of the smaller body is generally 
not constrained and therefore does not help with mass 
determination. Restricting the spin of the smaller mass 
does not help because at high mass ratios (and therefore 
NS-like m 2 ), the contribution to the PN phase from X 2 
is suppressed. The leading order spin corrections enter 
the PN phase with magnitude so X 2 ’s influence to 

the phase evolution is down-weighted relative to xi by 
O ((mass ratio)^). 

5. DISCUSSION 

In this paper we investigated the capability of the Ad¬ 
vanced LV network to distinguish between NSs and BHs 
from the inspiral-only waveforms using a large pop¬ 
ulation of plausible detections. This is the first large- 
scale study to characterize compact binary PDFs includ¬ 
ing spin-precession effects over a broad range of masses, 
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Fig. 6 . — Component mass recovery colored by the spin of 
the lower mass. Open circles mark the true values for nine 
representative binary systems from the simulated population. 
The spin of the lower mass is not well constrained, and does 
not strongly correlate with the mass parameters. Therefore, 
restricting the spin of the smaller body to be consistent with 
our prior expectations does not impact our ability to distin¬ 
guish between black holes and neutron stars. Both x- and 

y-axes are in units of Mq . „ .11 , r . 

mass ratios, and spins. (Jur study does not factor m 

systematic effects from real detector noise or differences 
between template waveforms and the true gravitational 
wave signal. Our study is limited to inspiral-only 
waveforms for simulation of signals and template 
waveforms for recovery because available precess- 
ing merger/ringdown models are not valid over 
the full prior volume. For many of our simulated 
signals the merger will be detectable and may 
help improve component mass estimates. Fur¬ 
ther improvements may come from including PN 
amplitude corrections. 

We arrive at four main conclusions from our analysis: 

1. When are we certain of at least one NS? For most 
systems with m 2 < 2 Mq, and all systems with m 2 < 
1.5 Mq. 

For larger-mass neutron stars (2 < m 2 < 3) the 90% 
credible intervals frequently extend into the low-mass BH 
regime. This tendency for high-mass neutron stars in 
NSBH systems to be consistent with low-mass black holes 
poses a challenge for determining the maximum NS mass 
from LV observations alone. 

2. When are we certain of a BBH with both masses above 
3 Mq ? When the mass of the smaller object exceeds 
~ 6 Mq. 

The m 2 posteriors for BBH signals seldom reach below 2 
Mq, so if a maximum neutron star mass were indepen¬ 
dently confirmed to be mmax ~ 2Mq, then TV’s ability to 
discriminate between BBH and BNS/NSBH populations 
would be significantly improved. 

3. When is a black hole definitely not in the mass gap? 
When component mass m 2 > 10 Mq the 90% credible 
intervals do not reach into the mass gap. 

It may prove challenging to confirm its existence because 
NSs with masses above ^ 2 Mq have error bars reach¬ 
ing into the m 2 € [3,5] Mq interval from below while 
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BBH systems with m 2 ^ 10 Mq can have credible in¬ 
tervals that extend below 5 Mq. Our findings suggest 
that inferring a bimodal mass distribution based solely 
on the observed sample of compact binaries especially 
without reliable predictions for the underlying compact 
object mass distribution will be challenging. 

4. When are we certain of a mass-gap object? Only in 
rare circumstances when the binary is nearly edge-on 
and spins are oriented in the orbital plane. For typical 
binaries component mass errors are larger than the 
mass gap. 

We simulated an additional 100 injections with both 
masses in the mass gap between [3,5] Mq. We found 
> 95% of the systems had some posterior support for 
a NSBH system. The majority of plausible mass gap 
sources yield credible intervals that are similar to, or ex¬ 
ceed the width of the mass gap itself (2 Mq). Only 
sources in low-probability alignments (spins in 
the orbital plane and/or edge-on orientations) 
were constrained to be 3 < m < 5 Mq. However, 
m 2 posteriors from mass gap BBHs do not reach below 2 
Mq . While we will not be able to say with any certainty 
that a given source occupies the mass gap, we will at least 
be alerted to the fact that the binary contains either an 
unusually high-mass neutron star, or a pair of unusually 
low mass black holes, with either possibility providing 
plenty of intrigue for the astrophysics community. 

Assuming flat mass distributions, we used simulated 
LV detections to infer the relative fraction of NSs and 
BHs in and above the mass gap, and estimate how 
many detections would be needed to confidently con- 


J. Abadie et al. (The LIGO Scientific Collaboration), (2014), 
arXiv: 1411.4547 [gr-qc] 

F. Acernese et al. (Virgo Collaboration), 

Class.Quant.Grav. 32, 024001 (2015), arXiv:1408.3978 [gr-qc] 

E. E. llanagan and 1’. Hinderer, Fhys.Rev. D77, 021502 (2008) 
arXiv:0709.1915 [astro-ph] 

J. S. Head, C. Markakis, M. Shibata, K. Uryu, J. D. Creighton, 
et al, Phys.Rev. D79, 124033 (2009), arXiv:0901.3258 [gr-qc] 

B. D. Lackey, K. Kyutoku, M. Shibata, P. H. Brady, and J. L. 
Friedman, Phys.Rev. D89, 043009 (2014) 

arXiv: 1303.6298 [gr-qc] 

L. Wade, J. U. E. Creighton, E. Ochsner, B. D. Lackey, B. E. 
Farr, T. B. Littenberg, and V. Raymond, 

Phys. Rev. D 89, 103012 (2014), arXiv:1402.5156 [gr-qc] 

T. Li, W. Uel Pozzo, S. Vitale, C. Van Den Broeck, M. Agathos, 
et al., Phys.Rev. D85, 082003 (2012), arXiv: 1110.0530 [gr-qc] 
L. Sampson, N. Cornish, and N. Yunes, 

Phys.Rev. D87, 102001 (2013), arXiv:1303.1185 [gr-qc] 

J. Abadie et al. (LlCU Scientihc), 

Class. Quantum Grav. 27, 173001 (2010) 
arXiv:1003.2480 [astro-ph] 

F. Ozel, D. Psaltis, R. Narayan, and J. E. McClintock, 
Astrophys.J. 725, 1918 (2010), arXiv: 1006.2834 [astro-ph.GA] 

W. M. Farr, N. Sravan, A. Cantrell, L. Kreidberg, C. D. Bailyn, 
et al., Astrophys.J. 741, 103 (2011) 
arXiv: 1011.1459 [astro-ph.CA] 

L. Kreidberg, C. D. Bailyn, W. M. Earr, and V. Kalogera, 
Astrophys.J. 757, 36 (2012), arXiv:1205.1805 [astro-ph.HE] 

K. Belczynski, C. Wiktorowicz, C. Fryer, D. Holz, and 
V. Kalogera, Astrophys.J. 757, 91 (2012) 
arXiv:1110.1635 jastro-ph.CA| 

C. Kochanek, Astrophys.J. 785, 28 (2014) 
arXiv:1308.0013 [astro-ph.HE] 


elude that gap BHs exist. Depending on whether the 
underlying mass distribution is heavily dominated by 
low- or high-mass systems, ten to hundreds of detections 
are needed to confirm objects in the gap. Assuming a 
plausible population-synthesis model where BNSs, NS- 
BHs, and BBHs are well separated in chirp-mass space, 
iMandel et al\ (|2015l) found that a mass gap was statisti¬ 
cally distinguishable with a only few tens of detectionsit 
is clear that forecasts for what we may learn from LV ob¬ 
servations of compact binaries are subject to large vari¬ 
ance while uncertainty about the true mass distribution 
persists. Continued effort in theoretically understanding 
that distribution, and how LIGO/Virgo observations can 
be used to test such theories, will be of great value as we 
begin assembling a catalog of compact binaries coales¬ 
cences. 

6. ACKNOWLEDGMENTS 

We thank Atul Adhikari, Claudeson Azurin, Brian 
Klein, Brandon Miller, Leah Perri, Ben Sandeen, Jeremy 
Vollen, Michael Zevin for help running the MCMC analy¬ 
sis. Will Farr, Carl-Johan Haster, and Ilya Mandel pro¬ 
vided important comments about our calculations and 
results. TBL and VK acknowledge NSF award PHY- 
1307020. BF was supported by the Enrico Fermi Insti¬ 
tute at the University of Chicago as a McCormick Fellow. 
SC thanks the US-UK Fulbright Commission for per¬ 
sonal financial support during this research period. DEH 
acknowledges NSE CAREER grant PHY-1151836, the 
Kavli Institute for Cosmological Physics at the Univer¬ 
sity of Chicago through NSE grant PHY-1125897, and an 
endowment from the Kavli Eoundation. Computational 
resources were provided by the Northwestern University 
Grail cluster through NSE MRI award PHY-1126812. 


D. Clausen, A. L. Piro, and C. D. Ott, 

Astrophys.J. 799, 190 (2015), arXiv: 1406.4869 [astro-ph.HE] 

I. Mandel, C.-J. Haster, M. Dominik, and K. Beiezynsk, 

(2015), 10.1093/mnrasl/slv054, 
arXiv:1503.03172 [astro-ph.HHJ 

L. Bianchet, Living Hev. Hei. 9, 4 (2006). 

P. C. Peters and J. Mathews, Phys. Rev. D 131, 435 (1963) 

B. Sathyaprakash and B. Schutz, Living Hev.Hei. 12, 2 (2009), 
arXiv:0903.0338 [gr-qc] 

C. Cutler and B. B. Blanagan, Phys. Hev. D 49, 2658 (1994), 
arXiv:gr-qc/9402014 

E. Poisson and C. M. Will, Phys.Rev. D52, 848 (1995) 
arXiv:gr-qc/9502040 [gr-qc] 

M. Hannam, U. A. Brown, S. Fairhurst, C. L. Fryer, and 1. W. 
Harry, Astrophys.J. 766, L14 (2013), arXiv:1301.5616 [gr-qc] 

K. Chatziioannou, N. Cornish, A. Klein, and N. Yunes, 
Astrophys.J. 798, L17 (2015), arXiv:1402.3581 [gr-qc] 

J. Veitch, V. Raymond, B. Barr, W. 4'arr, P. Cratt, et aL, 
Phys.Rev. D91, 042003 (2015), arXiv:1409.7215 [gr-qc] 

A. H. Nitz, A. Lundgren, D. A. Brown, B. Ochsner, U. Keppel, 
et al, Phys.Rev. D88, 124039 (2013), arXiv:1307.1757 [gr-qc] 

K. -W. Lo and L.-M. Lin, ApJ 728, 12 (2011), 
arXiv:1011.3563 [astro-ph.HBj 

A. Buonanno, B. Iyer, B. Ochsner, Y. Pan, and 
B. Sathyaprakash, Phys Rev D 80, 084043 (2009), 
arXiv:gr-qc/0907.0700 

M. Hannam, P. Schmidt, A. Bohe, L. Haegel, S. Husa, et al, 
Phys.Rev.Lett. 113, 151101 (2014) arXiv:1308.3271 [gr-qc] 

S. Vitale, R. Lynch, J. Veitch, V. Raymond, and R. Sturani, 
Phys.Rev.Lett. 112, 251101 (2014) arXiv:1403.0129 [gr-qc] 

J. Aasi et al (LIGO Scientihc Collaboration, Virgo 
Collaboration), (2013), arXiv:1304.0670 [gr-qc] 






7 


M. Vallisneri, Phys. Rev. D 77, 042001 (2008) 
arXiv:gr-qc/0703086 

I. Mandel,'Phys. Rev. D 81, 084029 (2010) 


